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We propose a theoretical model for branching instabilities in 2-dimensional fracture, offering 
predictions for when crack branching occurs, how multiple cracks develop, and what is the geometry 
of multiple branches. The model is based on equations of motion for crack tips which depend only 
on the time dependent stress intensity factors. The latter are obtained by invoking an approximate 
relation between static and dynamic stress intensity factors, together with an essentially exact 
calculation of the static ones. The results of this model are in good agreement with a sizeable 
quantity of experimental data. 



I. INTRODUCTION 

The phenomenon of crack division, i.e. the splitting 
of a single primary crack into two or more branches, 
whose dynamics develops independently, is studied in 
thin plates of different materials (glasses, plastics, metals 
etc.) pj. A classical example from more than thirty five 
years ago is shown in Fig. 1, which exhibits a crack 
pattern observed in a thin araldite plate 0- In this 
example the cracks go throughout the thickness of the 
sheet. This is to be distinguished from apparently sim- 
ilar side branching instabilities in which side branches 
appear in addition to the main crack, when the veloc- 
ity of propagation exceeds a critical value 0, ■ There 
exists an important difference between the two phenom- 
ena: in the former cases all the branches crack through 
the plates (some time referred to as "macrobranching" ) , 
whereas in the latter experiments, near the onset of the 
instability, the side branches have a width which is con- 
siderably smaller than the thickness of the plate (some 
time referred to as "microbranching" ) . Notwithstand- 
ing attempts to interpret the latter phenomenon using 
2-dimensional theories 0, |(| , it appears that the inter- 
action of crack fronts of different thickness necessitates 
a 3-dimensional theory which is daunting at present. In 
this paper we limit our discussion to a 2-dimcnsional the- 
ory that pertains only to plates with cracks going through 
the plate. 

The aim of this paper is to develop an approximate 
theory of crack dynamics, including crack bifurcations 
and multiple crack competition. More specifically, we 
address the following questions: 




FIG. 1: Crack pattern in an araldite tensile sheet 0. 



1. Given a straight propagating crack in a 2- 
dimensional material, when and how the first bi- 
furcation occurs? The bifurcation event itself was 
studied successfully in a recent paper by Adda- 
Bedia Q . In our approach we are able to examine 
the dynamics of the bifurcated cracks. 

2. Given a bifurcated crack, what is the stability of a 
symmetric branched configuration? We show that 
there exists an instability towards geometric per- 
turbations, making one branch growing on the ex- 
pense of the other which gets arrested. 

3. What is the geometry of multiple (two or more) 
branches? Obviously the interaction between mul- 
tiple branches results in curved cracks. The model 
that we propose is able to follow the dynamics and 
the resulting geometry. 

4. How consecutive bifurcations come about, and 
what are the resulting crack patterns? 

5. Can one apply the 2-dimensional theory devel- 
oped here to microbranching? The answer will 
be shown to be negative; there are crucial 3- 
dimensional aspects of microbranching that need 
the 3-dimensional theory in its full galore for an 
appropriate treatment. 

In order to answer all these questions, we take the point 
of view that the dynamics of the tip of each crack is 
determined by the elastic field in its very vicinity. Near 
each tip one expands the stress field as usual, 
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Here v is the instantaneous tip velocity, {r, 9} are polar 
coordinates at the crack tip and t is time. K 1 (t) and 
K u (t) are the dynamic stress intensity factors, and the 
functions S are known universal functions of 9 and v. 

The central element in our model is the adoption of the 
Hodgdon-Sethna equations for the crack tips Q. These 
equations were derived in the context of quasi-static crack 
propagation, and were shown to be in agreement with 
quasi-static experiments in Here we employ these 
equations in the dynamic context, invoking the results 
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and their comparison with experiments for justification. 
Consider a local coordinates system located at the crack 
tip, in which t and n denote the tangential and normal 
directions respectively. For the crack tip location r tlp we 
write the equations 
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Here v is the instantaneous crack tip velocity, and / is 
a positive material parameter. Needless to say, these 
equations appear simpler than they really are. The ac- 
tual calculation of the dynamic stress intensity factors for 
evolving cracks of complicated geometry is very far from 
trivial. The bulk of Sect. |H]is devoted to the presenta- 
tion of approximate schemes to compute these objects. 
The main idea is to relate the dynamic stress intensity 
factors to their static counterpart 0, and then to com- 
pute the latter using the method of iterated conformal 
maps that had been presented for fracture problems in 
recent papers 0>EI- This material is reviewed in Sect. 

nn 

Having at hand the stress intensity factors one can use 
them in Eqs. 10) for each crack tip. In Sect. IIVI we de- 
scribe how this is done for multiple branch dynamics. We 
demonstrate that the theory is successful in describing 
crack bifurcations, crack arrest, and successive bifurca- 
tions. The geometry can be studied in great detail and 
compared (very favorably) with available experiments. 
Finally, in Sect. |V|we consider the applicability of this 
theory to microbranching. As mentioned above, the con- 
clusion is negative. Sect. I VII offers a summary of the 
paper and some concluding remarks. 



II. DYNAMICS OF MULTIPLE CRACKS 

To implement Eqs. (J2J we need first to compute the 
dynamic stress intensity factors, and second the velocity 



v at the tip of each branch in a multiple branch configu- 
ration. We start with the stress intensity factors. 



Estimating the Dynamic Stress Intensity 
Factors 



The first task is the calculation of the dynamic stress 
intensity factors for a branched configuration. To this 
aim we invoke the formalism developed in 9], in which 
the dynamic stress intensity factors were related to their 
static counterparts. Based on some specific examples an 
approximate form of the dynamic stress intensity factors 
was obtained as a product of the static stress intensity 
factors and if " (of the instantaneous frozen configu- 
ration) and universal functions of the instantaneous ve- 
locity. For each mode of fracture one writes 



K u (t) 
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k u (v)K^(t) 



(3) 



We note that both the dynamic and the static stress in- 
tensity factors are considered time-dependent. For the 
static objects this dependence means freezing the actual 
configuration that is obtained at time t. The universal 
functions of the velocity are given by 



k x (v) = S 
k u (v) = S 

and S(Q is given by 
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Here cr, Cd and c s are the Rayleigh, dilatational and 
shear wave speeds respectively. 

This approximation was shown in the classical theory 
of fracture mechanics to be essentially exact for semi- 
infinite straight cracks |l2j. In our earlier work we ap- 
plied this approximate methodology for describing inter- 
acting large and small cracks [9j. The separable form 



of the dynamic stress intensity factors is trivially correct 
for very small velocities, since the functions ki and ku 
tend to unity for v — > 0. For finite velocities the sepa- 
rable form is not exact, but we expect it to yield good 
approximations when the velocities are small fractions 
of the typical wave speed and the typical distance be- 
tween the evolving tips is small. We expect this since the 
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mutual information about the location of each tip is car- 
ried by waves; under the specified conditions these waves 
can deliver the required mechanical information before a 
substantial change in the state of the system took place. 
We will show below that these conditions are satisfied 
in the early stages of the branching instability, therefore 
allowing us to use the separable form quite confidently 
in the present context. It is difficult to quantify a-priori 
the range of validity of the approximation. Therefore 
we invoke the final results and their agreement with var- 
ious experimental results to support the quality of the 
approximation. 

Clearly, this approximation is a huge simplification, 
calling for solving the static equilibrium field equations 
rather than the full dynamical field equations. Of course, 
one still has to face the difficult problem of static non- 
trivial geometries, but this problem was solved quite gen - 
erally using the method of iterated conformal maps |10| , 
and demonstrated in the context of complex crack ge- 
ometries in [llj. 



B. The velocity of the crack tips 

To close Eqs. as a consistent mathematical system 
we need to compute the velocity of each tip in terms of 
the dynamic stress intensity factors. The basic idea is to 
employ the energy balance that equates the energy re- 
lease rate into the crack tip region (denoted by G) to the 
dissipation involved in the crack propagation (denoted 
by r). The classical theory of linear elasticity fracture 
mechanics |T^ provides the energy release rate into each 
tip region: 



(6) 



where E and v are the Young's modulus and Poisson's ra- 
tio respectively and A 1 (v) and ^4 n (i>) are universal func- 
tions given by 



(1 - u) eg D(v) 

V 2 y/1-V 2 /C 2 s 



(7) 



An(v) c 2 D(v) ' 

D(v) = 4 ^(1 - v 2 /c 2 ) (I - v 2 /c 2 s ) - (2 - v 2 /c 2 s f . 

Note that D(v) vanishes at the Rayleigh wave speed, v — 
±cr. 

For concreteness, consider a two-dimensional infinite 
medium loaded at infinity with a uniform constant tensile 
stress a™. A long straight crack propagates at an instan- 
taneous velocity V; at some critical velocity, when the 
crack length is L, the crack bifurcates into two branches 
of lengths {£;} <C L with tip velocities defining 
angles {Ai7r} with respect to the direction of the crack 
prior to the bifurcation. Freezing the crack just at the 



bifurcation we denote its static stress intensity factor as 

K<f>\ 
i 

At each tip of the bifurcated crack we define the nor- 
malized stress intensity factors 
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Equating G to T we can rewrite Eq. @ at each branch 
tip as 
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Under the assumption that T is velocity independent we 
observe that the left hand side of Eq. © contains only 
material parameters. Therefore, a similar equation holds 
for the pure mode I crack propagating with velocity Vb 
just before the branching event, 



E r 

~ ~ 2 



K (o) 



(11) 



We conclude that the instantaneous velocity Vi of each 
crack tip is determined by Vb according to 



9 I (V b )=g I (v i )F?+g n (v i )F? I , Vz 



(12) 



Note that it would be better to index the stress intensity 
factor with an index i to stress that each tip contributes 
its own equation to the set. We avoid it in order not 
to overbear the notation. Notwithstanding, in case for 
which the dissipation function T becomes velocity depen- 
dent, we should introduce T(vi)/T(Vb) on the LHS of Eq. 
(JT2J. Bearing in mind that g,iyb) is a decreasing func- 
tion, the result of such a change would be a reduction in 

VbU- 

On physical grounds one seeks solution of these rela- 
tions for non- negative branch velocities Uj. In ref. it 
was shown that under the assumption that the branches 
start their evolution quasi-statically (i.e. vi — 0) a solu- 
tion appears first for symmetric branching at Vb = v c « 
0.475c s . For higher velocities Vb > v c one can have 
bifurcations in which the branches start off at a finite 
velocity. 

Finally, we bring the equations of motion to their final 
form as used below. Let be the angle between the 
tangential unit vector and the a;-axis. Denoting the tip 
position of the straight crack at bifurcation as rb, 1% = 
l r / P — r b\- We rewrite now the tangential and normal 
unit vectors at the tip of the crack in terms 9 and a re- 
scaled time tK^f — > t. These changes transform Eqs. 
into 



Mi 
dt 
d(k 
dt 



(13) 
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These equations, in conjunction with Eqs. i|12|) . define 
our dynamical system. Note that the velocities are mea- 
sured in units of c s . We selected v = 0.25 which implies 
Cd = Vic s and cr = 0.9194c s . It is important to notice 
that the set of equations for each branch tip is coupled 
to the equations for the other tips via the functions F T 
and F n . All the results presented below are obtained by 
the following procedure: for each instantaneous branched 
configuration the functions F 2 and F Jf are calculated us- 
ing the method of iterated conformal mappings, then the 
velocities of the tips are calculated using Eqs. (|12[1 and 
finally the increments in length and angle are calculated 
according to Eqs. (fH!|) . 



III. STATIC BRANCHED CONFIGURATIONS 

A crucial ingredient in our model is the calculation 
of the static stress intensity factors for an arbitrary 
branched crack configuration. The general approach to 
this problem, based on the method of iterated conformal 
mappings, is presented in all detail in ref. |lflj . The es- 
sential building block is the composition of the conformal 
map from the exterior of the unit circle to the exterior of 
the complicated crack shape, using a functional iteration 
of a fundamental conformal map that adds one single 
bump to the unit circle. This scheme enables us to solve 
for the entire stress field for any branched crack config- 
uration. The static stress intensity factors are extracted 
from the near tip fields using the method explained in 
[TfT | . In the present context the crack tips are dressed by 
a finite curvature determined by the size of the bump of 
the fundamental map which is added at each iteration. 
While appropriate for comparison with most realistic ex- 
periments, where some blunting of the crack tip is always 
present, this finite curvature means that comparison with 
mathematical models with infinitely sharp cracks should 
be done with care. 

In the literature there is only limited amount of works 
presenting calculations of stress intensity factors for in- 
plane problems with branched crack configurations. In 
ref. [5j such a calculation for an infinitesimal symmet- 
ric branched configuration is provided using a numerical 
solution of an integral equation. In order to ascertain 
the reliability and accuracy of our calculation we con- 
sider first a similar configuration, choosing the length I 
of each branch such that £/L = 0.5 x 10~ 3 . We could 
not select a smaller ratio due to the finite curvature of 
the crack tip (in addition, our numerical scheme which 
is based on truncated Fourier expansions cannot deal ef- 
ficiently with minute geometric details, since the series 
truncation becomes inaccurate). The normalized stress 
intensity factors F r and F n as a function of A, where 2irX 
is the angle between the branches (see inset) are shown 
in Fig. |2 This figure should be compared with Fig. 4 in 
Ref. [f| . It is clear that the figures are in good agreement 
(the locations of the maximum of F I and zero crossing 
of F zl are nearly identical), although there is a slight 
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FIG. 2: The normalized stress intensity factors F t and F TI as 
a function of A, where 2n\ is the angle between the branches 
(see inset). The ratio of branch length and the main crack is 
e/L = 0.5 x 10" 3 . 

overestimation in F t due to finite size effects. 

Second, we considered an asymmetric branched con- 
figuration in which both branches have the same length 
IJL = 0.5 x 10~ 3 , while one of them is located in the di- 
rection of the main crack and the other creates an angle 
of 7rA relative to that direction. The normalized stress 
intensity factors F t and F u for both tips as a function of 
A are shown in Fig. To our best knowledge there is no 
calculation available in the literature for this configura- 
tion. Since mode III models are usually in a qualitative 
agreement with their in-plane counterparts we present in 
the appendix the calculation for a mode III asymmetric 
branched configuration and present the resulting normal- 
ized stress intensity factors in Fig. Indeed, the mode 
I component of the in-plane calculation shows the same 
qualitative behavior as its mode III counterpart. 

IV. THE DYNAMICS OF MULTIPLE 
BRANCHES 

In Ref. H it was found that under the conditions 
that the branches start off quasi-statically (with zero 
velocity) and with K u w 0, the critical velocity is 
Vb = v c w 0.475c s and the branching angle A7r is 0.137T. 
Note that A is determined by the zero crossing of K u pre- 
sented in Fig. [21 The velocity is "critical" in the sense 
that it is the first velocity for which branching is ener- 
getically possible. Since there is no specification of the 
mechanism of branching, one should treat it as a lower 
bound for the branching velocity. For every Vb > v c one 
can find a solution with the same A (since it is deter- 
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FIG. 3: The normalized stress intensity factors F } and Fjj 
for both tips as a function of A, where 7rA is the angle between 
the branches (see inset). The ratio of branches length and the 
main crack is £/L = 0.5 x 1CP 3 . 
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FIG. 4: The mode III normalized stress intensity factor for an 
asymmetric branched configuration. 7rA is the angle between 
the branches and the ratio of the branches length and the 
main crack length is £/L — 0.5 x 10~ 3 . The upper (lower) 
curve corresponds to the forward direction (inclined) branch. 

mined by the independent condition K It ss 0), but with 
a non-vanishing velocity of the branches. Indeed, exper- 
imentally it appears that the branches do not emerge 
quasi-statically as implied by the solution of 5]. Note 
that as long as Vb is not much larger than v c , the branches 
velocities are relatively small. Bearing in mind that the 
distance between the tips is also relatively small, we can 
use the separable form of the dynamic stress intensity 
factors quite confidently. In this section we analyze the 
post-branching dynamics for various physical conditions. 



FIG. 5: The velocities of the branches as a function of time 
for Vb ~ 0.50 and Vb = 0.55. The resulting dynamics are such 
that the unperturbed branch is arrested, while the perturbed 
one returns after a short time to the original crack path. A 
representative resulting crack pattern is shown in the inset. 
See text for details. 



A. Stability Analysis 

Motivated by the experimental evidence that symmet- 
ric branches do not emerge quasi-statically, we study the 
stability of the symmetric configuration against small ge- 
ometrical perturbations. We consider the possibility that 
dynamical instabilities prevent the development of the 
branching event even though it is energetically allowed 
at Vb — v c . Consider a symmetric branched configura- 
tion with l/L = 0.5 x 1(T 3 and A = 0.13. Introduce 
a positive small perturbation 81 = 0.5 x 1Q~ 2 £ to the 
length of one of the branches, and integrate the dynam- 
ical Eqs. (|12I13(I for various Vb > v c . Note that as Vb 
increases so does the velocity of the emerging branches. 
Figure presents the resulting dynamics for Vb — 0.50 
and Vb = 0.55. A representative resulting crack pattern 
is shown in the inset; the unperturbed branch competes 
with the perturbed one and eventually dies out (i.e. it 
gets arrested due to screening effects). The velocities of 
the branches are plotted as a function of time. The time 
to arrest can be identified as the point where the veloc- 
ity of the unperturbed branch vanishes. By comparing 
the data for the two branching velocities, it is clear that 
the time to arrest increases substantially as Vb increases, 
therefore we deduce that the increment of instability de- 
creases with increasing Vb- We conclude that symmetric 
branched configurations are unstable against small geo- 
metrical perturbations at least for the branching veloci- 
ties we considered. We should stress that larger Vb's are 
not considered here due to the expected deterioration in 
the quality of the approximation embodied in the sep- 
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FIG. 6: The relative change in length A£/£ as a function of Vb 
for 51 = 2 x 10 _2 I (lower curve) and 81 = 0.5 x 10" 2 ^ (upper 
curve). Note that the continuation of the lines intersects the 
x-axis at the critical branching velocity Vb = v c « 0.475, 
below which branching is energetically forbidden. 



arable form of the dynamic stress intensity factors. In 
this regime it is reasonable to believe that other dynamic 
effects are important and might stabilize the symmet- 
ric configuration. The resulting dynamics are such that 
the unperturbed branch is arrested, while the perturbed 
one returns after a short time to the original crack path. 
Therefore, the main effect of these attempted branching 
events is the sudden deceleration of the crack. We suggest 
to interpret the instability as a possible explanation for 
the fact that in macrobranching events the branches do 
not emerge quasi-statically since then the configuration 
is very sensitive to perturbations and probably cannot 
be observed on macroscopic scales. On the other hand, 
for larger branching velocities, for which the branches 
emerge with finite velocities, both branches coexist, they 
are less unstable to small perturbations and therefore can 
grow to observable sizes. 

In the framework of stability analysis we also consider 
the final length of the arrested branch. Denoting by A£ 
the difference between the final branch length and its 
initial length £, we show in Fig. the relative change 
in length A£/£ as a function of Vb for two values of the 
fixed ratio 5£/£. The dependence seems to be approx- 
imately linear for both. Note that the continuation of 
the lines intersects the x-axis at the critical branching 
velocity Vb — v c ~ 0.475, below which branching is ener- 
getically forbidden. In passing, we note the resemblance 
of Fig. to the experimental results obtained in § for 
the microbranching instability. There it was found (see 
Fig. 3 in ||) that the branch length increases approxi- 
mately linearly with the mean crack velocity and vanishes 
for the critical one. 




FIG. 7: The crack pattern that was formed by three successive 
branching events. The branching velocity was set to Vb = 0.5, 
i.e. slightly more than v c . We introduced a positive small 
perturbation 5£ = 2 x W~ 2 £ to the length of the lower branch. 
Upon reaching the branching velocity again we introduced 
a symmetric configuration with the same £ and A with no 
perturbation added by hand. The bar shows the scale of the 
process relative to the initial crack length L. 



B. Successive Branching Events 

In light of the observation of asymmetric branch 
growth with one of them arrested, we follow now the 
evolution of the surviving branch. This branch then ac- 
celerates to the critical branching velocity and may bi- 
furcate again. In this subsection we study the patterns 
formed by such multiple successive branching events. In 
Fig. we present the crack pattern that was formed by 
three successive branching events. At the first event we 
introduced a positive small perturbation 5£ — 2 x 10~ 2 £ 
to the length of the lower branch. At the second event 
we introduced a symmetric configuration with the same £ 
and A with no perturbation added by hand. We stopped 
the evolution of the system at the onset of the third event 
since within our numerical precision we could not deter- 
mine its outcome; in some cases it turned out that the 
upper branch outruns the lower and in others vice versa. 
The results shown in Fig. are in qualitative agree- 
ment with the experimental results shown in Fig. |H1 We 
note here that the seemingly up-down anti-correlation be- 
tween the winning branches of successive events in Fig. 
13 is reminiscent of the spatial ordering observed in the 
microbranching instability (see Fig. 4c in This sim- 
ilarity, in conjunction with the resemblance discussed in 
relation with Fig. [BJ suggests that our 2-dimensional 
theory might have captured some of the features of the 
microbranching instability, although it cannot be directly 
applied in that case (see Sect. IvT) . 

The crack pattern shown in Fig. \7\ reveals a definite 
time scale r which separates successive branching events. 
We can understand this time scale with the following ar- 
gument: at branching, the length L of the crack is much 
larger than any other length scale like £ or 81. Accord- 
ingly, the stress intensity factors that determine the ve- 
locities are very weak functions of £. On the other hand, 
we expect the pattern of events to be of a linear size of 
the order of £. For a fixed S£ we thus expect r to be of 



7 




FIG. 8: A crack pattern in a tensile 260 x 120 x 5 mm araldite 
plate Q]. The plate was notched symmetrically at the edge. 
Due to minute asymmetries in the production of the notches, 
only one of them propagates. A second branching event can 
be observed after some time. In this event almost symmetric 
branches emerge from the branching point and coexist until 
one outruns the other and then curves towards the symmetry 
line. Attempted branching events can be clearly observed 
before the successful branching event took place. All these 
features are in pleasing correspondence with our discussion. 
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FIG. 9: Symmetric branching dynamics for different branch- 
ing length L with V& = 0.5. The plot shows the various crack 



patterns in rescaled coordinates. A power law y 
C, = 0.8 was added as a guide for the eye. 



with 



the order of I. Therefore, for a general perturbation 8£ 
we expect 

r~i,(!f) ■ 

In passing we note the relevance of the prediction l|14fl to 
experiments. In real materials there are always inhomo- 
geneities and asperities that introduce the perturbations 
that were modelled here by 81. Denoting the density of 
such asperities by p, we then estimate the length £ before 
meeting the first asperity by p . The size of the as- 
perity will determine 81. It appears worthwhile then to 
test the relation (|14fl in light of this identification. 

C. Symmetric Branching 

As discussed above, at higher branching velocities, 
symmetric branches can coexist for a longer time, and 
it is interesting to determine the typical profiles of such 
symmetric branches. Consider an experiment in which a 
crack of length L in a long strip of width W bifurcates 
and two symmetric branches of length £ emerge. As the 
branches start propagating with K u m 0, they cannot 
change their direction as long as £ -C L, W . On the other 
hand, when I grows to the order of the smaller between 
L and W, the branches will curve. In our theory there is 
only one length L (the system is infinite) and we study 
the curving of branches of length comparable to L. 

Figure shows several branching scenarios for different 
branching length L. It is tempting, after the example 
of Q, to fit power laws to these profiles. The different 
profiles can be approximated in the limited range that we 
consider by a power law y ~ with 0.7 < C < 0-8- The 



fit in the figure corresponds to £ = 0.8 and was added 
as a guide for the eye. It should be immediately said 
however that the profiles are neither universal nor true 
power laws. They represent a transient behavior between 
two straight lines. Initially the branches start off with an 
angle \ir = 0.137T. Finally there is an asymptotic fixed 
angle that depends on the geometry of the system. In an 
infinite medium this final angle satisfies < A < 0.13, 
while in a strip of finite width A = 0, i.e. the branches 
propagate eventually parallel to the boundaries of the 
strip. 

As an example of the relevance of our calculation to 
actual crack branching events, Fig. QJJI shows two shadow 
photographs of different stages of a symmetric branching 
event in a glass plate . The left panel shows the onset of 
branching with straight branches as long as £ is smaller 
than L and W. The right panel shows the developed 
branching configuration. The strong similarity with Fig. 
is obvious. 



V. HOW ABOUT MICROBRANCHING? 

With the relative success of the model proposed here 
in reproducing the dynamics and geometry of macro- 
branches it is of course tempting to see whether also 
the geometry of microbranches can be gleaned from the 
present 2-dimensional theory. To this aim we have at- 
tempted to follow the dynamics of a side branch in an 
asymmetric configuration in which the main branch is 
forced to emerge in the forward direction, see the inset 
of Fig. |3J We selected A = 0.2 such that K I1 of the side 
branch is approximately zero. This initial configuration 
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FIG. 10: Two shadow photographs of a symmetric branching 
event in a 300 x 100 x 9 mm glass plate p|. The left panel 
shows the early stages of the branching process where the 
branches are almost straight, while the right panel shows the 
developed curved branching configuration. 
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FIG. 11: A typical resulting configuration for asymmetric 
branching. The side branch is almost immediately arrested, 
while the main branch temporarily deflected from its straight 
path. The net effect of such an event is a strong fluctuation 
in the velocity of the main branch. The bar shows the scale 
of the process relative to the initial crack length L. 



was introduced as a constraint on the system to reflect 
the local symmetry breaking observed in experiments . 
Generally, this asymmetry results in a velocity difference 
between the side branch and the main branch. 

Figure 1111 shows a typical resulting configuration in 
which the velocity of the main branch fluctuates in a 
way similar to the experimentally observed velocity fluc- 
tuations, while the side branch was almost immediately 
arrested. Since we were mainly interested in the dynam- 
ics of the side branch, i.e. its trajectory and life time, 
we concluded that a model that treats the main macro- 
scopic branch and the microscopic side branch on equal 
footings is doomed to fail. The meaning of this result 
is that the energy flux into the microscopic side-branch 
is dramatically underestimated in a 2-dimensional model. 
We tried to increase artificially the energy release rate G 
into the near tip region of the side-branch by a constant 
factor, such as to increase its initial velocity. Neverthe- 
less, the velocity of the side branch dropped immediately 
to zero, reflecting the huge screening effect of the main 
branch. Therefore, although some of our results in Sect. 
IIVI show similarities with various features of the micro- 
branching instability, we propose that the phenomenon 
of microbranching is essentially 3-dimensional and can- 
not be modelled directly by a 2-dimensional theory. 



VI. SUMMARY AND CONCLUSIONS 

In summary, we have introduced dynamical equations 
of motion for crack tips which depend only on the stress 
intensity factors at the tips. Generally speaking, the cal- 
culation of these objects is daunting. By adopting an ap- 
proximate separable form for the dynamic stress intensity 
factors in terms of their static counterparts and universal 
velocity dependent functions we achieved a huge simpli- 
fication that results in tractable dynamics. Instead of 
complicated field equations we can reduce the theory to 
ordinary differential equations for the crack tips. Com- 
plex events like crack bifurcations, branch competition, 
branch arrest and successive bifurcations are studied in 
detail and compared with experiments. The good news is 
that the comparison with experiments is very encourag- 
ing. The bad news is that the quality of the approxima- 
tion cannot be easily assessed from first principles. One 
expects that for low velocities and small distances be- 
tween the multiple crack tips the approximation should 
be quite good. What is the range of validity can at this 
point be gleaned only from comparison with experiments, 
and these are relatively old and not detailed enough, 
maybe giving the false impression of a good agreement. 
It thus seems very worthwhile to conduct new experi- 
ments in light of this theory in order to test it in some 
detail. The gained transparency and simplicity of the 
theory seems a very good motivation for such an under- 
taking. 
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APPENDIX A: STATIC MODE III 
ASYMMETRIC BRANCHING 

The aim of this appendix is to derive the static stress 
intensity factors at the tips of an asymmetric branched 
mode III configuration. Since the mode III problem is 
described by Laplace equation for the displacement field 
u z (x,y) 

V 2 u z (x,y)=0, (Al) 

the solution of this problem is readily given if the confor- 
mal map, z = $(w), from the exterior of the unit circle to 
the exterior of a branched crack configuration is known. 
The general formalism for obtaining such a map is known 
in the literature for a long time [2j. Here we adapt the 
general formalism to the problem at hand and solve it. 
Consider the following map 

- Aj~ 1 («;-e 4oi ) Al (to-e io,a ) As («;-e <a8 ) A » , (A2) 
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where A is a real constant, < ot\ < u-i < < 2n 
and A 3 = 2 — Ax — A 2 . The points {e tak } are mapped to 
the origin and therefore are branch points. For a^-i < 
9 < a k the phase of <5(u/) is fixed by the crack branch 
angle and a local maximum of |$(w)| is obtained at e %l3k . 
The parameters of the map, i.e. A, {otk} and {/3fe}, can 
be calculated by demanding that for w = 1 arg(z) = 
and by the conditions $'(e l/3fc ) = 0, |$(e i/3fc )| = l k . 
Here {£ k } are the lengths of the crack branches. These 
three conditions can be translated into the following set 
of equations 



E oijXj = 2tt 



A (0Lj- fik 

Xj cot 1 



= 



4A n 



sm 



/5fe - Q.j 



(A3) 



which can be solved numerically. Having the conformal 
map for the required configuration at hand we are mainly 
interested in the stress intensity factors for the various 
crack tips. The solution for Eq. (|A1| is given by 



2 f i. 



ip(z) + if>{z) 



(A4) 



where ip(z) is an analytic function. The stress compo- 
nents in polar coordinates are given by 



a rz (r,9) - ia ez {r,9) = e l \'[re w ) = 



(A5) 



Here the tilde denotes the usual transplantation. In the 
near vicinity of a given crack tip z k we use the following 
expansion 

a rz (r,9) = -§2L sin (9/2) + 0{y/Pj 
V27rr 



a 8z {r,9) = -^cos(9/2)+0(V^) ■ (A6) 
V 2nr 

These two equations can be rewritten as 

-ie ie ' 2 K, 



a rz {r,9) - i(Te z {r,9) = 



Si + 0{VP) ■ (A7) 



On the other hand, near the crack tips we can expand 
the conformal map to obtain 

z-z k = $(w) - $(u k ) ^ i$"(w fc )(^ ~ ^fe) 2 



$'(w) ~ - w fc ) = V2*"(w fc )(z - z fc ) 

(A8) 

The last expression, in the light of Eq. I|A5(I , shows the 
explicit relation between the square root singularity of 
the stress field near the crack tip and the derivative of 
the conformal map. Let us denote z — z k = re l i s k^+s) 
with 5 k = Ylj—i Xj and consider the direction tangent to 
the crack tip, i.e. 9 = 0. By comparing Eq. (|A5|) with 
Eq. (|AT|) we obtain 



K 1U = i<p'(wk) 



$"(w fe )e i,5 >= 7 



(A9) 



With this result at hand we can calculate the stress in- 
tensity factor at each tip since the solution in the w-plane 
is known to be 



-icr 



(A10) 



where cr^ is the applied stress at infinity. Now we are in 
a position to analyze infinitesimal asymmetric branched 
configurations. We choose t\ = L, £ 2 — £3 — £, Ai = 1 — A 
and A2 = A with £/L = 0.5 x 10 -3 . The resulting stress 
intensity factors are presented in Fig. 0|m the text. 
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